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Abstract 

We propose a model which explains the experimental observation of spon- 
taneous spin polarization of conducting electrons in quasi-one-dimensional 
AlGaAs/GaAs channels H],^. We show that a ferromagnetic order is a 
generic property of a quasi-one- dimensional conducting channel embedded 
in a Wigner crystal. We also discuss gate voltage, magnetic field and temper- 
ature dependences of the channel's conductance. 
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Recently both theories and experiments on conductances of quasi-one dimensional con- 
ducting channels attracted a lot of attention (See for example Quasi-one-dimensional 
confining potential has been made by applying split gates to two-dimensional electron gas 
in AlGaAs/GaAs heterostructures (See the insert in Fig.l). In this case electron motion is 
quantized in x direction. If the confining potential changes adiabatically in "y" direction, 
we can neglect backscatterings of electrons incident to the channel from the two-dimensional 
reservoirs. For noninteracting electrons at zero temperature and in the absence of an external 
magnetic field H, the conductance of the system is quantized 



Here ki = 0, 1,2... is an integer associated with number of quantized levels in x direction, 
which are filled with electrons. The prefactor 2 in Eq.l reflects the spin degeneracy of 
electron states a.t H = 0. In experiments the electron concentration in two-dimensional 
semiconductor structures was controlled by a gate voltage Vg. Eq.l would manifest itself in 
the existence of plateaus in the V^-dependence of GiVg) quantized in units When the 
concentration of electrons in the channel is high, or A;i ^ 1, Eq.l is in very good agreement 
with experimental data ll|-§]. However, at small concentration of electrons, when ki ~ 1, 
experimental data [§,0 deviate from Eq.l dramatically. The following anomalies have been 
observed: 

a). At small temperature T, in addition to plateaus at G = ^,2^, the dependence 
G{Vg) exhibits "smeared plateaus" at G ~ 0.7^ and G ~ 

Furthermore, in the noninteracting electron picture, an external magnetic field H parallel 
to the film lifts the spin degeneracy and each plateau splits into two plateaus quantized in 

2 

units ^. Again, this is in agreement with experiments at large ki. However, the plateaus at 
G = 0.7^, ^1.7^, 2^ do not split in the external magnetic field |^,^. It has been noted 
that the H dependence of positions of these plateaus can be explained if one assumes 
a ferromagnetic electron ground state in the channel at i/ = 0. 




(1) 



b). The first and third plateau 
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and fourth plateau (G = 2^). 

c) The plateaus at G = ^,2^ become weaker as T increases, while the plateaus at 
G ~ 0.7^, 1-7^ become stronger. 

Existing theories of one-dimensional interacting electron gas predict no spin polarization 
in the ground state and can not explain the measured dependence G{Vg, H,T) when the 
concentration of electrons in the channel is small. 

In this article we propose a qualitative explanation of experimentally observed anomalous 
features. Consider the confining potential V{x, Vg) = V'o(Fg)(l + -j^jrj) as shown in Fig.l. 
If the potential V{x) is shallow (the value of A is big) and the concentration of electrons in 
the channel is small, the ground state of the electron system is a Wigner crystal. In Fig.l 
positions of locahzed electrons in the Wigner crystal are indicated by sohd dots. We will 
assume that the Wigner crystal is pinned and the conductance of the system is determined 
by hoppings of point defects in the crystal (interstitial electrons or vacancies) , concentration 
of which is exponentially small. The activation energy E'^(x) of interstitial electrons is 
nonuniform and has a minimum £'^(Vo) = E^[x = 0) in the middle of the channel. 

At a critical value Vq — Vqc, melting of Wigner crystal occurs near x — and a quasi-one 
dimensional conducting channel is formed. Below we consider two scenarios where creation 
of the channel corresponds to a first and second order phase transition. In both scenarios 
the channel is ferromagnetic. 

I. SECOND ORDER LIFSHITZ PHASE TRANSITION 

Let us assume that at the critical value Vq = Vcq, for the first time, the activation energy 
of the interstitials becomes zero in the middle of the channel, i.e. El^{Vo — Vqc) — 0. This 
can be considered as a second order Lifshitz phase transition. The interstitial electrons are 
indicated in Fig. 1 by the symbol " * " . This transition is a precursor of melting of a crystal 
into a liquid. The existence of zero point motion defects in the middle of the channel is a 
nonuniform analogy of the model of supersolid proposed by A.F.Andreev and I.M Lifshitz in 
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the context of quantum melting of helium crystals 0. Such a state has never been observed 
in bulk helium crystals. Latter it was noted |0 that conditions for existence of such a state 
are less restrictive at a crystal-liquid surface, where the amplitude of quantum fluctuations 
of the lattice sites is larger. This leads to a quantum rough state of the surface, where 
zero point motion steps on the surface exist in the ground state, which, in turn, leads to 
the existence of surface crystallization waves at T = 0. Crystallization waves at small 
T, have been observed [§], while existence of solid-liquid helium quantum rough surfaces in 
d=3 case is still an open question P,[To||. In two dimensional case the existence of quantum 



rough surfaces is more plausible than in the three dimensional case because the amplitude 
of quantum fluctuations of particles at the surface is larger. 

To begin, we first consider the case when all spins are polarized in | direction, for 
example, by an external magnetic field parallel to the film. In this case the interstitial's 
energy spectrum in the quasi-one dimensional channel has a form 

^l{K,Py) = ^o{K + ^) + A| cos(apy). (2) 

Here = 1,2.. is an integer corresponding to quantization of the wave function in x 
direction, py is a momentum in y direction, eo = 2m* > ^ ~ ^2^/^^^)^ ^ characteristic 
dimension of the wave function in x direction, m* ~ -i^-^ is the effective mass of the 
interstitial, Aj^ is the band width for interstitials with spin J, direction, and a = ^ and 
n are interelectron distance and electron concentration in the Wigner crystal respectively. 
Band structures for interstitials moving in a polarized Wigner crystal in the case eo > Aj 
are shown in Fig. 2. 

In the absence of an external magnetic field H = 0, the interstitial electrons wave 
functions and spectrum inside the conducting channel depend crucially on the magnetic 
structure of the Wigner crystal. The magnetic structure of a uniform Wigner crystal is 
not known (See a discussion in We will assume that the magnetic exchange energy 

between electrons J has an antiferromagnetic sign and the ground state corresponds to 
either an antiferromagnetic state or a resonating valence bonds (RVB) state |T^. We will 



show that in both cases the conducting channel should be ferromagnetic. It has been noted 
|1^,[13| that in some range of concentrations the ferromagnetic state of the Wigner crystal is 
also possible. In this case, our results are valid at relatively high temperatures T > J. We 
believe that in experiments [^,^ this inequality took place because in the Wigner crystal 
J~ ImK 0. 

It has been known for a long time ||l5|-p0[l that a point defect embedded in a uniform 



antiferromagnetic quantum crystal is surrounded by a ferromagnetic region with a radius 

R ^ a{^y/^ » a. (3) 

We assume that Aj ^ J. The reason for the existence of the magnetic polaron is that A| 
is much larger than the energy band width of the defect embedded in an antiferromagnetic 
crystal. A competition between negative delocalization energy of the point defect inside the 
ferromagnetic polaron ( of order of 2_r^L* ) positive spin polarization energy of the 

magnetic polaron, (of order of ^(-f )^) leads to Eq.3. 
If T > J one should substitute J with T in Eq.3. 

Thus there are two characteristic lengths in the problem, L and R. Let us consider the 
case L < R when the magnetic polarons have sizes R and L my and x direction respectively. 
If N, the one-dimensional concentration of the interstitial electrons in the channel, is small 
(A^ <^ R~^), the magnetic polarons do not overlap; the magnetization M{y) is nonuniform; 
the channel conducts poorly and the conductance of the channel increases gradually with 
N (or with Vg). At iV > R-^, the magnetic polarons overlap and form a conducting 
magnetically polarized channel. The width of the magnetically polarized region in this case 
is of order of L. At this point the concentration of interstitials in the channel is small and, 
due to exchange spin splitting, they have spins |. The conductance of the channel exhibits 
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the first plateau at G = ^. 

The situation arising at a high concentration of interstitials N depends on a ratio between 
eo and Aj^. If eo ^ Aj, further increasing of N will lead to filling of the band k = 1, |. 
Since A| < A^, the kinetic energy of interstitials with spins f is small (See Fig. 2). Here 



A| is the width of the band with k = The width of the magnetic stripe does not 

increase as the band A; = 1 , t is being filled and interstitials in the band = 1 , t move in 
the uniform exchange potential associated with the uniform magnetization M inside the 
channel. Therefore the conductance exhibits the second plateau at G = which is not 
smeared. 

At last, let us discuss the temperature dependence of the first plateau, which corresponds 
to the band A; = 1, |. Suppose that (Vq — Vqc) < and |Vo — Voc\ is small and, consequently, 
the interstitial activation energy E'^ > is also small. In this situation interstitials thermally 
activated in the channel increase the magnetic moment M and make the distribution of M{z) 
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more uniform. As a result, the first plateau at G ~ ^ becomes stronger as temperature 
increases. Since i?^ is small, one can neglect thermal fluctuations of the Wigner crystal 
lattice. This is different from solid ^He case where the activation energy of point defects is 
large. Therefore the crystal melts before point defects contribute significantly to magnetic 
properties of the sample. 

The presented above model explains the following experimental facts. 

a) . It predicts the ferromagnetic ground state of the conducting channel and the existence 
of an additional step in the V^-dependence GiVg, H=0) at G = ^. 

b) .It predicts that the first plateau in the V^-dependence of G{Vg) is more smeared than 
the second one. 

c) .It explains why the first plateau at G ~ ^ becomes stronger as T increases. 

d) .It also follows from the presented above picture that at small concentration of electrons 
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in the channel when G ^ ^, and at small T, the system should exhibit giant negative 
magnetoresistance in a parallel magnetic field. This is because the radius of the magnetic 
polaron R{H) increases with H. Obviously, this effect takes place also in the variable range 



hopping regime [21 



At last, let us mention difficulties of the presented above model. 

1). Experimental data exhibit the plateau at 0.7^. It is not clear to us whether it 
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is possible to attribute the difference between 0.7^ and ^ to smearing of the first plateau. 

2). There are experimental indications that the first and the third plateau disappear 
at T = [^,^]. To this end, we would like to mention that the filling of the first 
band, as Vg changes, is a complicated process which involves overlapping of the magnetic 
polarons. Therefore, even in the approximation when there is no direct interaction between 
the interstitials, there is an indirect long range many-body effective interaction between the 
interstitials, which is mediated by spin degrees of freedom of the Wigner crystal. As a result, 
if is fixed, the spatial phase separation of the interstitials takes place at relatively small 



[jl8[. Recently this phenomenon has been reconsidered in the context of HTC The 
apparent differences between [^] and the situation discussed here are: a)the interstitials are 
neutral objects because they are screened by Wigner crystal lattice and b) in experiments it 
is the chemical potential rather than that is fixed. The latter means that at T = there is 
no phase separation. Instead, N{Vg) dependence exhibits a jump from zero to a finite value 
of the interstitial concentration = Ac > 0. The final state of the conducting channel in 
this case is uniform. As it is usual for first order phase transitions, the states with = and 
A = Ac are separated by an energy barrier. Thus, the formation of the conducting channel 
is a collective phenomenon which involves tunneling of many electrons into the channel. 
Therefore, at this point, we can not rule out a possibility that the G(yg,T, H) dependence 
is hysteretic and that the experimental results depend on the history of the temperature, 
the magnetic field and the gate voltage. 

3). In the approximation of noninteracting interstitials and at < cq the band = 1, f 
will be completely filled before the the chemical potential reaches the top of the band k = 1, |. 

2 2 

This would lead to reentrance from the plateau at G = 2^ to the plateau at G = ^ as 
increases. 

It is worth emphasising, however, that the presented above consideration is of a single 



We are grateful to D.Khmelnitskii for turning our attention to this difficulty. 
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particle nature. It neglects interstitial-interstitial and interstitial- Wigner crystal lattice in- 
teractions and requires detailed knowledge of interstitials band structures at relatively high 
interstitial energies. In experimental situation, = {^Y^'^ is not such a big number and 
the interaction effects could be very important at large A^. In particular, the quasiparticle's 
bandwidth could depend on A^. 

Some of these difficulties may be resolved if we consider the opposite limit eg ^ A^. 
In this case k* ~ levels with spin | should be filled before the chemical potential of 

2 

interstitials reaches the band k = 1,1. Therefore the conductance is quantized in unit ^ 
and plateaus in G{Vg) follow in a proper sequence. It is also plausible that the first plateau 
becomes stronger as T increases. It is not clear, however, why in the framework of this 
picture the third plateau at G = ^ should be more smeared than the second one at 

In principle, it would be possible to distinguish eo > and eo < A^ cases by different 
Vg{H) dependences. Here Vg are the gate voltages at which the chemical potential reaches 
the bottom of the i^^ band and a new plateau occurs. In the case Aj^ < eo, the gate voltages 
Vg and Vg correspond to bottoms of /c = 1, | and A; = 1, f bands. Therefore they move 
in the opposite directions as H increases. In the case A^ > eo, gate voltages Vg and Vg 
correspond to bottoms of k = 1 [ and k = 2 I bands and they move in the same direction 
as H increases. In experiments , Vg and Vg move in the same direction. It is not clear, 
however, whether this can also be explained as an orbital magnetic field effect which exists 
due to a finite width of the two dimensional electron gas 0,^. 

II. CONDUCTING CHANNELS WITH LARGE WIDTH 

In a uniform situation, the two-dimensional Wigner crystal melting is believed to be 
a first order phase transition, which occurs at a critical electron concentration ric. In a 
nonuniform case shown in Fig.l, a classical picture would predict that the transition takes 
place at a larger concentration of electrons in the middle of the channel n(x = 0) = > Uc 
because of existence of the crystal-liquid surface energy. The width of the conducting channel 
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at the point of transition is = (^)^/^ > a. Here 7 = n /^-^^^"|^~^^^"^^ ^^(n) and 
are energy densities of 2D Fermi liquid and Wigner crystal respectively, and a is the 
surface energy per unit length. 

Now consider a situation when the width D is large and the concentration of electrons 
in the channel is high. In this case we can view the system as a stripe of Fermi liquid 
sandwiched between two antiferromagnetic (or RVB) Wigner crystal (See Fig.Sa). 

It has been pointed out in the context of solid-liquid ^He interface that the ferro- 



magnetic ordering is a generic property of a quantum rough surface. It originates from the 
competition between kinetic energy of defects on the surface and the antiferromagnetic ex- 
change energy in the bulk of crystal. Since it is not known whether the Wigner crystal- Fermi 
liquid surface is in a quantum smooth or in a quantum rough state, we will consider the case 
of a quantum smooth boundary, which is less favorable for the ferromagnetism, and present 
a model which exhibits boundary ferromagnetism even in this case. 

We will model the system with the help of the Hubbard Hamiltonian defined on a semi- 
infinite two-dimensional square lattice shown in Fig.Sa. 



Hd = tY^ aj^ttja + c.c. + UY^ aj^aiial^ai^. (4) 



Here i,j labels adjacent sites on a square lattice, Uiaja]^ are annihilation and creation 
electron operators on a site i, a is a spin index labeling | and J, directions of electron spins, 
t and U are parameters which describe intersite tunneling and intrasite electron repulsion 
respectively. In the case of half-filling and at t ^ [/, the Hamiltonian describes a Mott 
insulator, which is an analog of the Wigner crystal. The spin ground state of the system is 



antiferromagnetic (or RVB) with exchange energy of order J ~ ^ 



We will describe the semi-infinite two-dimensional Fermi-liquid with the help of the 
Hamiltonian for noninteracting electrons = J2q ^gbqbq, where bq, 6^ are annihilation and 
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creation electron operators in the liquid, = ^ ~ Ep is the electron energy, q and m 
are the momentum and mass respectively, Ep is the liquid's Fermi energy. The tunneling 
between the Mott insulator and the Fermi liquid is described by the Hamiltonian 
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HT^tiJ2aibl + cc. (5) 

i',q 

Here index i' labels the localized states in the insulator at the boundary. The total Hamil- 
tonian is H — Hd + + Ht- In the case of half filling, the chemical potential of the Mott 
insulator hcs in the middle of the gap between the valence and the conduction band (See 
Fig. 3b). For a Wigner crystal, it means that the energy of creation of a vacancy and an inter- 
stitial electron would be the same. In general, there is no reason to expect such a symmetry 
in a real Wigner crystal. Therefore we introduce into the model a small concentration of 
impurity levels, which pin the position of the chemical potential below (or above) the middle 
of the gap. These impurities are marked by symbols "+" in Fig. 3b. The distance between 
the impurities will be much larger than all other characteristic lengths in the problem. To 
be concrete we consider the situation when the chemical potential is pinned near the bottom 
of the valence band and the energy for transferring an electron into the liquid and creation 
of a vacancy in the Mott insulator \sUi <tiU (See Fig. 3b). 

The tunneling between the liquid and the Mott insulator leads to a negative correction 
to the ground state energy of the system. The value of the correction depends on the spin 
configuration in the Mott insulator. In the second order perturbation theory with respect 
to ti, we have an expression for the correction to the energy per unit surface length 

8E'^^ y (6) 

where a is the lattice spacing, e{5}(n) is the energy of a vacancy moving in a spin config- 
uration {S} in Mott insulator, n labels the vacancy eigenfunctions in the Mott insulator 
and Ly is the sample length in y direction. Suppose that in the Mott insulator near the 
surface there is a ferromagnetic region of the width Lj. (See Fig. 3a). Then we can estimate 
e{5}(n) ~ ^cosipyO) + k1 ■ H^re A = 2t is the vacancy energy bandwidth in a com- 
pletely spin polarized Mott insulator, is an integer corresponding to quantized motion 
of the vacancy in x-direction, py is the quasi-momentum of the vacancy in "^/" direction 
(n = kx,Py)- We assume that in the case of the antiferromagnetic ordering the bandwidth 
of the vacancy is small and it is almost localized. 
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Expanding Eq.6 with respect to |^ ^ 1 and taking into account the energy of the 
spin polarization in the stripe of the width L^, we get an estimate for spin configuration 
dependent part of the ground state energy 



SEfoi t]f_ JL^ 



Eq.7 has a minimum at 

Lx = ci\ 



(8) 



Therefore the system is magnetically polarized at the surface provided jjr^ ^ 1- 

The statement about existence of the magnetic polarization near the surface can be 
proved in a more rigorous way starting with usual perturbation theory 

AE*^ = 1 E [de^{-^rU^,^'). (9) 



Here labels the lattice sites at the boundary, T{i,i') labels paths of a vacancy on the 
lattice intersecting the boundary at sites i and z', which preserve an initial spin configuration 
of the system ( An example of such a path is shown in Fig.Sa), Ji{i,i') is a function of the 
distance between sites i and i', which reflects the Friedel oscillations in the Fermi liquid, np 
is the number of steps a vacancy hops along a path F. While obtaining Eq.9 from Eqs.4,5, 
we neglected all terms which involve double electron occupancies of lattice sites. The factor 
C(nr) depends on the spin configuration in the system. In the ferromagnetic case, due to 
Fermi statistics of electrons, it is equal to (—1)"'". In the limit ppa ^ 1, Ji = 6i^i' and we 
only need to take into account paths which enter and leave the Mott insulator through the 



same lattice site i = i'. At this point the problem is reduced to that considered in ||T5|,|16 



where it was proven that the ground state of the system with one vacancy and infinite U is 
ferromagnetic. A resummation of these ferromagnetic paths leads to Eq.6. 

If the surface is in a quantum rough state the ferromagnetic ordering at the surface is 
even more favorable energetically p5 |. 
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The presented above picture is relevant when the width D of the Fermi hquid stripe is 
large, which corresponds to large ki in Eq.l. Strictly speaking, in this case, each plateau 

2 

with given ki also splits into two plateaus quantized in unit ^. However, the amplitude of 
this splitting of plateaus should be small. 

Let us consider what happens as concentration of electrons in the channel and, con- 
sequently, D decreases. It has been mentioned above that the mean field theory, which 
neglects quantum fluctuations of the position of the crystal-Fermi liquid boundary, predicts 
the creation of the conducting channel as a first order phase transition. The question arises 
whether quantum fluctuations of electron positions in the Wigner crystal can make this 
transition to be a second order one. The amplitude of quantum fluctuations of the position 
of the crystal-liquid boundary depends on whether the surface is in a quantum rough or a 
quantum smooth state. If it is in a quantum rough state, we can describe deviations of the 
boundary from its equilibrium classical position by a variable X{y,t) and use the following 
effective action 

r rlX rlX 

S = J dydm^f + a{^Y + V{X)] (10) 

to estimate the amplitude of quantum fluctuations of the position of the boundary. Here 
V ~ ^X"^ corresponds to the energy difference between the solid and liquid phase as the 
boundary is moved by distance X, ^ = 7^- The second term in Eq.lO corresponds to the 
surface tension and the first term corresponds to acceleration of steps along the surface. 

We would like to mention that the situation described by Eq.lO is different from that 
considered in 0. The electron density profile n{x) in the channel is determined by the 
long range Coulomb force. Therefore, in our case there is no redistribution of the electron 
density in the liquid associated with oscillations X{y,t). Using Eq.lO we get an estimate < 
{6Xy) >~ ^J—^ log ^. Here the brackets <> stand for the quantum mechanical averaging. 
If < {dXy) >^ D^, then the one-dimensional melting can occur as a second order phase 
transition and we arrive at the picture considered in the previous section. We believe that, 
in principle, this inequality can take place because the transition occurs at ~ 39 where 
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at short distances the crystal differs from the hquid very httle. 

In conclusion we would like to make the following remarks. At finite temperature the 
ferromagnetic ordering of the quasi-one dimensional stripe will be destroyed by thermal 
fluctuations. Strictly speaking, at T = 0, the presented above arguments only prove existence 
of strong ferromagnetic correlations at "small" distances. At this stage we can not rule out 
a possibility that at "large" distances the ferromagnetism will be destroyed by quantum 
fluctuations of spin configuration in a 2-D Wigner crystal. Quantization of the conductance 
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of the channel in units ^ , however, persists as long as the length of the channel is shorter 
than the characteristic distance of the destruction of the ferromagnetism. 



It has been shown that the ground state of strictly one-dimensional systems can 
not be ferromagnetic. However, in the experiment the system is not strictly one- 
dimensional. We are unaware of any theorem which prevents existence of ferromagnetism in 
quasi-one-dimensional system of the type shown in Fig. 3. We leave this questions for future 
investigation. 

We are grateful to C. Markus, D. Thouless, S. Kivelson, D. Khmelnitskii, E. Lieb and M. 
den Nijs for useful discussions. Work of B.S. is supported by Division of Material Sciences, 
U.S. National Science Foundation under Contract No.DMR-9625370. F.Z. is supported by 
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FIG. 1. The qualitative picture of the confining potential V{x). Solid dots correspond to 
positions of electrons in the Wigner crystal. The symbol " * " corresponds to an interstitial electron. 
The insert shows the geometry of the conducting channel. 
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FIG. 2. A schematical picture of band structures of interstitial electrons in the case eo 3> 
The dashed line indicates the chemical potential when the channel becomes conducting. 
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FIG. 3. a. The geometry of the Mott insulator (MI)-Fermi Hquid (FL)-Mott insulator struc- 
ture. The Mott insulator and the Fermi liquid are indicated by symbols "MI" and "FL" respec- 
tively, r represents a path corresponding to a term in the perturbation sequence in Eq.9. b. 
Energy diagram for the Mott insulator-Fermi liquid boundary. 
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